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Abstract 

We study atomic chains with nonlinear nearest neighbour interactions and prove the ex- 
istence of fronts (heteroclinic travelling waves with constant asymptotic states). Generalizing 
recent results of Herrmann and Rademacher we allow for non-convex interaction potentials and 
find fronts with non-monotone profile. These fronts minimize an action integral and can only 
exists if the asymptotic states fulfil the macroscopic constraints and if the interaction potential 
satisfies a geometric graph condition. Finally, we illustrate our findings by numerical simula- 
tions. 
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1 Introduction 

Nonlinear Hamiltonian lattices like chains of interacting atoms or coupled oscillators are ubiquitous 
in mathematics, physics, and material sciences. The most famous example, and most elementary 
model for a crystal, is a chain of identical atoms that interact by nearest neighbour forces. In 
reminiscence of the pioneering paper by Fermi, Pasta, and Ulam |FPU55j one usually refers to such 
systems as FPU or FPU-type chains. 

Although FPU chains are quite simple lattice models they exhibit a rich and complicate dy- 
namical behaviour, and we still lack a complete understanding of their dynamical properties. A 
mayor topic in the analysis of FPU chains is therefore the investigation of coherent structures such 
as travelling waves and breathers. Travelling waves are highly symmetric, exact solutions to the 
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underlying lattice equation. They can be regarded as the fundamental modes of nonlinear wave 
propagation and provide much insight into the energy transport in discrete media. In this paper 
we aim in contributing to the general theory by studying fronts, i.e., heteroclinic travelling waves 
that connect two different constant states. 

The dynamics of FPU chains is governed by the lattice equation 

Xj = $'(xj +1 - Xj) - &(Xj - Xj-l). (1) 

Here Xj = Xj(t) denotes the position of the j th atom at time t, <3? is the interaction potential, and 
the atomic mass is normalized to 1. Introducing the atomic distances rj = Xj+i — Xj and velocities 
Vj = Xj we can reformulate flTJ) as 

rj = v j +i-v j , v j = &(r j )-&(r j -i). (2) 

A travelling wave is a special solution to ([2]) that satisfies the ansatz 

Tj{t) = R(j - at), Vj (t) = V(j - at). (3) 

R and V are the profile functions for distances and velocities, a denotes the wave speed, and 
tp = j — at is the phase variable. In dependence of the properties of R and V travelling waves 
come in different types. Wave trains have periodic profiles and are investigated in [FV991 IPP001 
IDHM 06]. They describe oscillatory solutions to JT]) and provide the building blocks for Whitham's 
modulation theory. Another important class of travelling waves are solitons (or solitary waves), 
where R and V are localized over a constant background state. The existence of solitons in lattices 
is a nontrivial problem and has been studied intensively during the last 20 years. We refer to 
[FW941 ISW97I IFM021 IPanOSl ISZ071 IHer09| for variational methods, and to [hSoOOl HJ05] for an 
approach via spatial dynamics and centre manifold reduction. 

In this paper we study fronts which have heteroclinic shape and satisfy 

lim R((p) = r±, lim V((p) = v± (4) 

with (r_, V-) 7^ (r + , v+). Fronts have attracted much less interest than solitons, maybe because 
they only exist if the asymptotic states satisfy some very restrictive conditions. In particular, $' 
must have at least one turning point between r_ and r+, and this excludes for instance the famous 
Toda potential. Nonetheless, fronts in FPU chains appear naturally in atomistic Riemann problems, 
see |HR10| for numerical simulations, and are important in the context of phase transitions. 

The first rigorous result about fronts we are aware of is the bifurcation criterion from [IooOOj . 
It implies that fronts with small jumps between the asymptotic states exist only if <3?' has a convex- 
concave turning point. Recently, the existence of fronts was proven by variational methods in 
|HR09j . The existence theorem therein does not require the asymptotic states to be close to each 
other but is restricted to convex potentials The proof relies on a Lagrangian action integral for 
fronts with prescribed asymptotic states and uses the direct approach to establish the existence 
of minimizers. A similar approach is used in |KZ09al IKZ09bj to prove the existence of fronts for 
sine- Gordon chains. 

In this paper we generalize the method from |HR09| and prove the existence of fronts without 
convexity assumption on Our main result can be summarized as follows. 

Theorem 1. Action minimizing front solutions to ([2]) exist under the following hypotheses: 

(i) The asymptotic states and the front speed satisfy the macroscopic constraints, which take the 
form of three independent jump conditions. 
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(ii) The potential satisfies the graph condition with respect to the asymptotic states. 
(Hi) Some technical assumptions are also satisfied. 

Moreover, there is no front without (i), and no action minimizing front without (ii). 

The assumptions in Theorem [1] will be specified below. The macroscopic constraints, see Lemma 
[21 are algebraic relations and link fronts to energy conserving shocks of the p-system, which is the 
naive continuum limit of FPU chains. In particular, they determine the wave speed a and imply 
that the asymptotic strains r_ and r + cannot be chosen independently of each other. The graph 
condition reformulates the area condition from [HR09] and requires that the graph of $ is below 
the shock parabola associated with the asymptotic states. Both the macroscopic constraints and 
the graph condition appear naturally in our variational existence proof and guarantee that the 
action integral is well-defined and bounded from below. 

Closely related to fronts are heteroclinic waves with oscillatory tails. These are travelling wave 
solutions to (|2|) which approach two different periodic waves for ip — > ±00. Such oscillatory fronts 
are used to describe martensitic phase transitions and to derive kinetic relations in solids [BCSOla, 
IBCSOlbl IAP071 rVailO] . The only available existence results, however, concern piecewise quadratic 
potentials, which allow for simplifying the travelling wave equation by means of Fourier transform, 
see |TV05l ISCC051 ISZ09] . It remains a challenging problem for future research to give alternative, 
maybe variational, existence proofs that cover more general chains. 

The paper is organized as follows. In £J2] we discuss the macroscopic constraints and normalize the 
asymptotic states. Moreover, we reformulate the front equation as an eigenvalue problem for a 
nonlinear integral operator. In we set the existence problem into a variational framework and 
characterize fronts as minimizers of an action integral. Or main technical result is Theorem 1161 and 
guarantees that this action integral attains its minimum on a suitable set of candidates for fronts. 
The proof uses separations of phases, which are introduced in £ )3.4I and allow to extract convergent 
subsequences from action minimizing sequences. Finally, we present some numerical simulations in 

m 

2 Preliminaries about fronts 

Substituting the travelling wave ansatz ([3]) into ([2]) yields 

v-fcR(<p) + + 1) - V(<p) = 0, v-fcV(<p) + &{R(<p)) - &(R(<p - 1)) = 0, (5) 

which is a nonlinear system of advance-delay-differential equations. Moreover, combining both 
equations we readily verify the energy law 

^ (^ 2 (V) + HR(<P))) + *'{R{<p))V{<p + 1) - &(R{<p - l))V(<p) = 0. (6) 
2.1 Macroscopic constraints for the asymptotic states 

We now derive the macroscopic constraints that couple the front speed a to the asymptotic states 
(r±, v±) from @. To this end we consider continuous observables ip = ip(r, v) and denote by 

hP(r, v)]\ := ip(r+, v + ) - ^(r_, vJ) and (tp(r, v)) := |(t/?(r_, u_) + ip(r + , v + )), 

the jump and mean value, respectively. 

The following result was proven in |HR09j (see also [AP07] ) by integrating ([5]) and (|6|) over a 
finite interval [-N, N] and passing to the limit N — > 00. 
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Lemma 2. The asymptotic states of each front satisfy 

a\[r}\ + |H| = 0, a\[v]\ + [*'(r)] = 0, a\[\v 2 + $(r)| + |[<D'(r>| = 0. (7) 

Heuristically, Lemma [2] reflects that fronts transform into shock waves when passing to large 
spatial and temporal scales. The jump conditions ([7]) precisely mean that the asymptotic states 
correspond to an energy conserving shock for the p-system and imply that each front satisfies 
mass, momentum, and energy. The p-system is the naive continuum limit of FPU chains under 
the hyperbolic scaling and reads 

d T r = d y v, d T v = d y &(r), (8) 

where r = et and y = ej denote the macroscopic time and space, respectively, and e > is a small 
scaling parameter. The conservation laws in ([8]) correspond to mass and momentum, and imply 
the conservation of energy for smooth solutions, that is 

d T (\v 2 + $(r)) = d T (v&(r)). (9) 

The jump conditions for ([9]), however, is independent of the jump conditions for ([8|). More details 
about the p-system and energy conserving shocks can be found in |HR10l lHR09j . 

Using the discrete Leibniz rule [V'lV^I = l^iK^} + (V-'iHV^I we readily verify that (J7|) implies 

[*(r)] = MK^Cr)), a 2 = W(r)]\/\[r]\. (10) 

Conversely, for any (r_, r + ) with (|l(J|) i there exist - up to Galilean transformations - exactly two 
solutions to (J7|) which differ in sgno". We now characterize the geometric meaning of (|10p and refer 
to Figure [1] for an illustration. 




Figure 1: To each front there exists a parabola that touches the graph of $ in both r_ and r + . 
Consequently, the signed area between the graph of and the secant connecting r_ and r + vanishes 
the stripe [r_, r + ]. 

Lemma 3. The following conditions are equivalent: 
(i) (a, r_, r+) fulfils 

(ii) there exists a parabola that touches the graph of & in both r_ and r + , 

(Hi) the signed area between the graph of $' and the secant connecting r_ to r + sums up to zero 
in [r_, r + ]. 

Moreover, each condition implies that has at least one turning point between r_ and r + . 
Proof. Consider the parabola f(r) = \ar 2 + br + c. The touching conditions 

f(r±) = *(r±), f'(r±) = $'(r±) 
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are equivalent to 



M^l+frM = [*(»■)], h(r 2 ) + b(r) + c=mr)), \[r]\a = \[&(r)]\, a(r) + 6 = <<I>'(r)>, 
and by ^|[r 2 | = (r)|[r]| we conclude that (i) and (ii) are equivalent via 

a = a 2 , b={&(r))-a 2 (r), c =<$(r)) - (r) + a 2 «r) 2 - I(r 2 )) . 

The equivalence of (ii) and (iii) is immediate since the secant has slope a 2 , and $' must have a 
turning point because otherwise the graph of <3?' would be either below or above the secant. □ 

Condition (|1U|) i is the kinetic relation for fronts and reveals that the asymptotic states cannot 
be chosen arbitrarily. More precisely, for given r_ and r + we can choose a and \[v ]\ such that the first 
two jump conditions in (which correspond to mass and momentum) are satisfied. However, for 
the energy condition ([7J)3 to hold, r_ and r + must additionally fulfil pOjb . Form this we conclude 
that fronts do not exist if is either convex or concave, and that in general we cannot prescribe 
both r_ and r + . 

We emphasize that ([7]) is in general not sufficient for the existence of fronts, i.e., there exist energy 
conserving shocks in the p-system that can not be realized by a front in FPU. In fact, it was proven 
in [IooOOj that fronts bifurcate from convex-concave but not from concave- convex turning points 
of This disproves the existence of subsonic fronts with small jump heights although there exist 
the corresponding energy conserving shocks. 

In order to prove the existence of action minimizing fronts we shall additionally to ([7]) require 
that the graph of $ is below the parabola defined by the asymptotic states, see Assumption In 
particular, our existence result provides a front for Example A from Figure [1] but does not cover 
Example B, see Remark 1121 and the examples in £jH 

2.2 Normalization and reformulation 

For our analysis in Sj3] it is convenient to normalize the asymptotic states and to reformulate the 
front equation ([5]) as an eigenvalue problem for a nonlinear integral operator. 

Lemma 4. Up to affine transformations we can assume that 

(j = l, r± = ±l, u± = =Fl, *'(±1) = ±1, $(±l) = i. (11) 

Moreover, with (|11|) the front equation is equivalent to 

l 

W = A&(AW), (AW)(<p)= J W(<p)d<p, (12) 

i 

where W is a normalized profile with lim (/ ,_ > -|- 00 W((f) = ±1. 
Proof. Let U and W be two normalized profiles such that 

R(V) = (r) + i|[r]|C% + l/2), V(<p) = (v) + \[v\W{<p). 
Using the first two jump conditions from ([7J) we readily verify that ([5]) transforms into 
-&U{<p) = W(cp + 1/2) - W(x> - 1/2), -tW{fp) = &(U(<p + 1/2)) - %'(U(<p - 1/2)), (13) 



5 



where the normalized potential 



$(u) 



s 1 ir 11 \ 2 ($'( r )) 1 4<$(r)) 
$((r) + i[r] u) - ±- > . / u + - - ,.J . L /. 



|[$'(r)]||[r]| V ' ^.l , , |[$/(r)]| 2 mr)][r] 



satisfies $'(±1) = ±1. Moreover, we have $(—1) = $(+1) = i if and only if the third jump 
condition 03 is satisfied. Towards (fl~2|) now suppose (fTT]l . Integrating (fT3|) i we find U = AW, 
where the constant of integration vanishes due to £7(±1) = W(±l) = ±1, and similarly we derive 
W = A$'(U) from (|T3j)9. □ 

The front parabola for normalized data (|lip is r i— > ^r 2 and each solution to (|12|) can be viewed 
as a perturbation of the shock profile 

+ 1 for (p < 0, 

W sh (^) = sgn<^ = <( for = 0, (14) 
-1 for if > 0. 

Notice that the residual of W s lu that is W s h — ^4$'(^4W s h), has compact support. 



We proceed with some preliminary remarks about the action of a front. Heuristically, the action 
density in the normalized setting is given by 

\W 2 - §{AW) = \W 2 - \{AW) 2 + ^{AW) 

with 

*(r) = \r 2 (15) 

so the action integral formally reads 

C(W) = J \W 2 - \{AW) 2 + V(AW)d<p. (16) 

R 

Notice that is just the difference between the front parabola and see Figure [IJ and that C is 
well defined as long as W approaches its asymptotic states sufficiently fast. A further possibility 
for defining the action integral was introduced in [HR09J for monotone W and relies on the relative 
action integral 

£( W ) = J {\ w2 - HAW)) - (±Wl - $(AW sh )) dip. 

R 

Both approaches are linked by C(W) = C{W) — C(W s h) and the symmetry of A, compare Lemma 
[7J formally implies 

dC(W) = dC(W) = W - A$'{AW). 

In £}3]we give a slightly different definition of C, see (|18p and (|23p . and establish the existence of 
minimizers. 



3 Existence of fronts 

In this section we assume that the asymptotic states and the potential are normalized by (jlip and 
show that the fixed point equation (|12p has a solution in some appropriate function space. 
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3.1 Assumptions 

We rely on the following standing assumptions on the function ^ from (|15p . Examples and coun- 
terexamples are given in £JH 

Assumption 5. ^ is continuously differentiable and satisfies the following conditions: 
(G) graph condition: \P(tt) > *(±1) = for all uel, 
(X) genericity: *"(±1) > and ^{u) > for all u ^ ±1, 

(M) monotone asymptotic behaviour: ^(u) is decreasing for u <C — 1 and increasing for m> 1. 

Condition (G) has a natural interpretation in terms of $ and can easily be reformulated for 
non-normalized data: It precisely means that the front parabola touches the graph of in both r_ 
and r + but is above this graph in all other points. Moreover, (G) is equivalent to the area condition 
from HH()9j , which characterizes the signed area between the graph of and the secant connecting 
r_ to r+. With positive and negative sign for above and below the graph of respectively, the 
area condition reads as follows. The signed area is non-negative in each stripe [r_, r] with r > r_ 
and non-positive in each stripe [r, r+] with r < r+. We refer to Figure [1] for illustration, where 
positive and negative area are displayed in dark and light grey colour, respectively, and recall that 
the signed area vanishes in the stripe [r_, r + ] due to (fTUlh . 

We mention that (G) is truly necessary for the existence of action minimizing fronts, see Remark 
[T2l The conditions (M) and (X), however, are made for convenience and might be weakened for 
the price of more technical effort. 

Remark 6. (X) is equivalent to 

(S) supersonic front speed: $"(±1) < 1, 

and (M) implies 

(I) invariant set for <£': There exits a constant T > 1 such that <£' maps [— T, V] into itself. 

Proof. (S) follows from the definition of ^ in f)15|) . Towards (I) we exploit (M) to choose T > 1 such 
that > u for u < — T and <3?'(u) < u for u > T. Then we set T = max{r, max^^p |^>'(ti)|}. □ 

3.2 Functionals and operators 

We denote by L p , W 1 '^ and C k the usual function spaces on the real line, abbreviate the L p -norm 
by || • ||p, and write 

{Wi, W 2 ) = [ W 1 {^)W 2 {ip)d^ 
Jr 

for the dual pairing of W\ € VP and W 2 € L p ' with 1 = 1/p+l/p'. 
Lemma 7. The averaging operator A has the following properties: 

1. A maps L p into L°° n W 1,p C C for all 1 < p < oo with 

(AW)'(<p) = W{<p+l)-W(<p-±) 
and \\AW\\ LP < \\W\\ LP , \\AW\\loo < \\W\\ LP , \\AW\\ w i, P < 3||W||up- 

2. A is symmetric in the sense that (AW\, W 2 ) = {Wi, AW 2 ) holds for all W\ G L p and 

W 2 e lp'. 
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3. A is self-adjoint in L 2 with spectrum spea4 = {g(k) : k € M} where g(k) = § sin (|) . 

Proof. The first two statements are straight forward. The third one follows since A diagonalizes 
in Fourier space via Ae lhlp = g(k)e lkip . □ 

We now introduce the afhne space 

n = {w : W-W sh el 2 }, 

where W s h is the shock profile from (|14|) . Exploiting Lemma the Taylor expansion of <£' around 
±1, and the properties of W s h we then find 

AW, &(AW), A&(AW) en, W - AW, W - A 2 W e L 2 . (17) 

for all W € T~L. In view of the action integral (|16p we also define a functional A4 on L 2 by 

M(^) = i j V 2 - (AV) 2 dp = i yV - ^ 2 ^)^ <ty>, 

R K 

and a functional A/" on % by 

jV(W) =M{W -W sh ) + \ j W 2 h - (AW sh ) 2 dp + J (W - W sh ) {W sh - A 2 W sh ) dep. (18) 

R R 

Notice that N(W) is well defined on H as both W 2 h — (AW^) 2 and W s ^ — A 2 W & \ 1 have compact 
support. Moreover, if W — W s h decays sufficiently fast for <p — > ±oo (say W - W sh e L 1 ), then we 
have 

M{W) = \ Jw 2 - (AW) 2 dp=\ j (W- A 2 W)Wd^. (19) 

R R 

Lemma 8. T/ie functional A4 is non-negative and weakly lower semi- continuous on L 2 . 
Proof. Denoting the Fourier transform of V by V we find 

M{V)= [ [l- g{k) 2 )V{k) 2 dk= \Wl-g 2 V\\ 2 



with g as in Lemma [7J This gives the desired result as V n — ^ V^, implies V„ — and hence 

Lemma 9. The functional M is Gateaux differentiable on T~L with derivative 

dJ\f(W) = W — A 2 W € L 2 . (20) 
Moreover, N is invariant under shifts in ip-direction, and satisfies 

M{W 2 ) = JV(Wi) + M(W 2 - Wi) + (W 2 - Wi, Wi - A 2 W X ) (21) 
for all W x , W 2 € U. 
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Proof. A direct computation with W G H and SW G L 2 shows 

(dM(W), SW) = {(W - W sh , SW) - {AW - AW sh , ASW)) + <^ sh - A 2 W sh , SW) 
= (W- W sh , SW) - (A 2 W - A 2 W sh , SW) + (W sh - A 2 W sh , SW) 
= (W-A 2 , SW), 

and this gives (|20p . Towards the shift invariance we approximate W by 

W n = X[ -n,n](W-W sh ) + W sh 

where X[—n,n] ^ s the indicator function of the interval [— n, n]. Then we use (|19p for W n to find 
£{W n ) = C(W n (- + ip)) for all shifts (p, and passing to the limit n — > oo gives the desired result. 
Finally, by definition we have 

N(W 2 ) - JV(Wi) = A^(W 2 - W sh ) - M(Wi - W sh ) + <W 2 - Wi, W sh - A 2 W sh ) 

and 

M(W 2 - W sh ) -M{W l - W sh ) = M{W 2 - Wx) + (W 2 - W u W x - W sh - A 2 W X + A 2 W sh ), 
so ()2ip follows from adding both identities. □ 
To conclude this section we consider the functional 

V(W) = J ^(AW)df, 



which gives the non-quadratic part of the action integral (|16p 

Lemma 10. V is well defined on H with 

c\\AW - sgn(AW) \\ 2 < V(W) < c\\AW - sgn(AW) || 2 (22) 

for some constants c and c that depend only on HAWHoo- Moreover, V is Gateaux differentiate on 
H with derivative dV(W) = A 2 W - A$'(AW). 

Proof. Let W G H be given and recall that U = AW satisfies U G % n L°° due to Lemma [71 
Condition (X) provides two constants c and c such that 

c(u — sgnn) 2 < ^(u) < c(u — sgnn) 2 

holds for all \u\ < HC/Hoo, and we conclude that V is well defined and satisfies (|22p . Finally, (|17p 
provides A 2 W — A<k'(AW) G L 2 , so both the existence of and the formula for dV follow from a 
direct calculation. □ 

3.3 Variational setting 

We now introduce the action functional on T~L by 

C(W) = M(W) + V(W). (23) 

In virtue of Lemma [9] and Lemma [TUI the functional C is well defined, shift invariant, and Gateaux 
differentiable with derivative 

dC(W) = A 2 W - A&(AW) G L 2 , 



9 



and we conclude that each minimizer of C in T~L must solve the front equation ()12|) . However, 
proving the existence of minimizers in % turns out to be difficult and therefore we restrict C to the 
convex subset 

c = {w e Hnw 1 -™ : ||w||oo < r, ll^'Hoo < 2r} . 

Notice that the ansatz W G C is reasonable due to condition (J) and since the front equation (|12p 
combined with Lemma [7] implies W G H n W 1,0 °. 

In order to link fronts to minimizers of £ in C we observe that the properties of $' and „4 guarantee 
C to be invariant under the L 2 -gradient flow of C. To see this we consider the explicit Euler scheme 

W ' y T\(W) = W - XdC(W) = (1 - \)W + \A$'{AW) (24) 

with small step size A. 

Lemma 11. The set C is invariant under the action of T\ for < A < 1. Consequently, each 
minimizer of C in C solves the front equation (|12|) . 

Proof. For W G C let P = <§>'(AW) and recall that AW, P, AP G H according to (HZ]). Combining 
(I) with \\AWWoc < T and (AP)' = P(- + 1/2) + P(- + 1/2) gives 

||P||oo, pp||oo<r, ||(XP)'|| QO <2r, 

and hence AP G C. Since C is convex we also have T\(W) G C for all < A < 1, and passing to the 
limit A —¥ we then establish the invariance of C under the L 2 -gradient flow of C. In particular, 
each minimizer of C in C must be a stationary point for the gradient flow of C and hence a solution 
to the front equation. □ 

To complete the existence proof for fronts it remains to show that C attains its minimum in C. 
We prove this in the next section by using the direct approach, that means we construct minimizers 
as limits of minimizing sequences. 

A particular problem we have to overcome in the subsequent analysis is that C is not coercive 
on %. In fact, as illustrated in Figure [2] there exist sequences (W n ) n C C with extending plateaus 
at —1 or +1. These plateaus contribute neither to M nor V but may imply 

\\W n - WshOn + Olla 00 

for all choices of the relative shifts ip n . Heuristically it is clear that the cartoon from Figure [2] 
cannot be prototypical for action minimizing sequences, but in order to proof this we need a better 
understanding of sequences with bounded action. 

We conclude with a remark about the necessity of the graph condition (G) and refer to ^5] for 
numerical examples. 

Remark 12. Suppose that ^ satisfies (M), (X), and \I/(±1) = 0, but violates (G) because there is 
some ii* with |it*| < T and \£(u*) < 0. Then C is unbounded from below. 

Proof. We define a sequence iW n ) n C % of piecewise linear profiles by W n (tp) = it* for \ip\ < n and 
W n (tp) = sgnip for \tp\ > n + 1/r. By construction, W n — A 2 W n is supported in I n U (—I n ) with 
In = [n — 1 , n + 1 + r] , and a direct calculation shows 

C{W n ) <C + 2(n- -oo, 

where C is some constant independent of n. □ 
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Figure 2: Sketch of a sequence with bounded action and two extending plateaus at ±1: graphs 
of W n and U n = AW n in Black in Gray, respectively. The shaded areas indicate corresponding 
separations of phases with I n j = ip n j + Jj as in Lemma [TBI 

3.4 Separation of phases for sequences with bounded V 

To characterize the qualitative properties of a profile U = AW with W € C we interpret U < \ 
and U > o as negative phase and positive phase, respectively, and regard intervals in which U 
takes intermediate values as transition layers. Obviously, adjacent plateaus of different height are 
separated by a transition layer and each U must exhibit at least one transition layer as it connects 
-1 to +1. 

We next exploit the uniform L°°-bound for W € C to derive a lower bound for the ^-contribution 
of each transition layer. To this end we introduce 

Zu = {<p: \U{(p)\<\}. 

which is nonempty, closed and bounded as U S C is continuous with U((p) — >• ±1 as <p — > ±oo. 

Remark 13. There exist constants f) > and p, > such that 

&(U(ip)) dip > p, 



<p-v 

for each U € C and (p G Z\j . 

Proof. This follows since U G C implies \U((p2) — U(ipi)\ < J Vl \U'(ip)\ dip < 2T\if2 — <pi\. In 
particular, we have | £/(</?) | < 3/4 for all — < fj = l/(8r) and the claim follows with ji = a/(2fj) 
and a = sup| u | <3 / 4 ^(u) > 0. □ 

In order to show that each function U possesses a finite number of transition layers we introduce 
the following definition. A separation of phases for a given profile U € C is a finite collection of 
closed intervals (transition layers) ii, I m , m > 1, such that 

1. the intervals are disjoint and ordered, i.e., 

min/i < maxii < min/2 < ... < max I m -i < min/ m < max/ m , 

2. Zjj is contained in I\ U ... U I m . 
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3. for each interval Ij there exists p>j G Z\j such that with [ipj — rj, pj + rj\ C Ij, 

Lemma 14. Xei W £ C be given and set U = AW . Then there exists a separation of phases 
I\, I m for U with 

m < floor(7 ? (W)//2) and 2fj < \Ij\ < 4P(W)fj/p, 

for all j = l...m. 

Proof. We define a finite number of points <p~j G Z\j and intervals Ij = (<pj — 2f), ipj — 2f\) iteratively 
as follows: p>\ is the smallest element of Zjj, i.e. U{p) < —\ for all ip < (p\. If Z\j \ I\ is empty 
we stop the iteration; otherwise we choose <p~2 to be the smallest zero of Z\j outside of I\. Then we 
have <p~2 — fti > 2f], so Remark 1131 yields 

V(W)> J V(U(<p))d<p+ J ^{U(p))dp>2fi. 

If possible, we now define (p^ as the minimum of Z\j \ [l\ U I2) and proceed iteratively until the 
iteration stops after m < floor (V(W)/fi) steps. By construction, we have Zy C UjLi^j an( ^ 
< 4f/ for all j but the intervals Ij may overlap. Finally, we obtain the desired separation of 
phases by merging overlapping intervals. □ 

Our main result in this section concerns sequences (W n ) n G C with bounded V(W n ). It guarantees, 
roughly speaking, the existence of compatible transitions layers which (i) have the same length, (ii) 
separate the same phases, and (Hi) depart from each other. For an illustration we refer to Figure 

12 

Lemma 15. Let (W n ) n C C be a sequence with limsup n _ s>00 V(W n ) < 00. Then there exists a 
not relabelled subsequence along with a finite number of intervals Ji, J m with the following 
properties: 

1. Each interval Jj is centred around zero, i.e., Jj = [—rjj, r/j] for some < rjj < 00. 

2. For each n there exist shifts <p n \ < ... < <p n ,m such that 

(a) p n ,i + J\i •••) (fn,m + J m is a separation of phases for U n = AW n , 

(b) Pn,j+i — fn,j — > 00 as n — )• 00 for all j = l...m — 1. 

3. There exists a choice of signs so, s m with sq = —1, s m = l, and Sj G {—1, +1} such that 

sgn(U n (tp)) = s for p < (p n ,i - m 

sgn(U n (p)) = si for p n ,i + m < p < p n ,2 ~ m 



n,m 



Sgn(U n (<p)) = S m -i for ip n ,m-l + Vm-l < p < <p. 
sgn([/ n M) = for (p ntm +r) m < pi 

hold for all n. 

Proof. Thanks to Lemma [H we can extract a subsequence such that each U n has a separation of 
phases that consists of m intervals I n> i, ...,I n ^ m where m is independent of n. We denote the centre 
of I n j by p n ,j, set J n j = I n j — p> n j, and notice that c < \J n ,j\ < c and p n ,j+i — Pnj > c for some 
constants c > c > independent of n and j. In particular, the intervals Jj = (J n J n j have finite 
length c < I Jj I < c. 

Our strategy for the proof is to refine to, the subsequence (W n ) n , the phase shifts p> n ,j, arid the 
intervals Jj in several steps. To this end we start the following algorithm at level 1. 
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Level k 



If m = k, then we stop the algorithm. 

If (f n ,k+i ~ Wfc — ^ oo as n — y oo along a subsequence, then we extract this subsequence 
and jump to level k + 1. 

If m < k and < sup n {Lp n ^+\ — tp n ,k) < oo, then we merge Jk and Jk+i as follows: At 
first we choose Jk sufficiently large such that 

Jk 3 Jk U ((fn,k+l - ¥n,k + Jk+l) 

n 

Secondly we define (p n ^k = p n ,k, rh = m — \ and 

Jj = Jj and if n j = (p n j for j < k, 

Jj = J j+ i and (p n j = (p n ,j+i for j > k. 

Finally we restart level k with ui f ip n j, and Jj instead of m, Lp n j, and Jj. 

This algorithm stops after a finite number of steps when m — k = 0. It provides intervals Jj and 
phase shifts (p n j for j = l...m and n £ N with lim n _ s . 0O v? n ,j+i — fn,j = oo for all j. By extracting 
subsequences we can also ensure that, for each n, the intervals (p nj + J n j are pairwise disjoint 
and provide therefore a separation of phases for U n . Finally, by extracting further subsequences if 
necessary we guarantee the existence of a choice of signs. □ 

3.5 Existence of minimizers for L 

We now finish the existence proof for fronts. 

Theorem 16. C attains its minimum on C and each minimizer is a front. 

Proof. Step 0. We start with some notations. For a given minimizing sequence (W n ) C C we 
define 

U n = AW n , S n = sgnlln, 

and for each K > we introduce the operator 

E K :C^C, E K W = X [-k,k](W - W sh ) + W sh . 

Here X[-k,k] denotes the usual indicator function, so we have H-E/fW — W\\2 — > as K — > oo for 
each W £ C. Finally, within this proof C always denotes a positive constant that is independent 
of n and K, but the value of C may change from line to line. 
Step 1. By assumption and M.{W n — W s h) > we have 

V(W n ) = C{W n ) - M(W n ) <C+\\W- PFshHocll^sh - ^ 2 ^ sh ||i < C. (25) 

Therefore we can extract (a not relabelled) subsequence for which Lemma [15] provides a finite 
number of intervals J\...J rru sequences of phase shifts (^> n ,j) n and a choice of signs (s n j) n . There 
exits at least one 1 < < m such that Sj„_i = —1 and sj t = +1, and since C is invariant under 
shifts we can assume that (f n ,j* = 0. With Jj f = [— r/j„, ?7jJ and due to lim n ^. 0O v? n ,j+i — ¥?n,j = oo 
we then have 

limsup U n (ip) < — 2 for tp < —T)j t , liminf U n {<p) > \ for ip > r]j t . 

n— ¥oo n—s>oo 
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By compactness we can extract a further subsequence such that W n — Woo weakly* in W In 
particular, W n converges to Woo uniformly on each compact interval, and hence 

E K W n E K Woo, E K U n EkUoo strongly in L 2 for all K. (26) 

Step 2. Towards Woo G C we show that W n —S n is uniformly bounded in L 2 . The first observation 
is that (|25p combined with (|22p implies 

\\U n -S n \\ 2 <C. (27) 
The second observation is that both S n — AS n and S n — A 2 S n are supported in the 1-neighbourhood 

Of In = [JT=A J j + 

<Pn,j)- Therefore, \I n \ < Ylj=i l^jl — C yields 
\\S n — AS n \\2 < C, 115,1 — .A 2 5 n ||i < C, |jV(5 n )| = \(S n — A 2 S n , 5 n )| < C, 
and by ([27]) we find 

||Z7 n — -4,5 n ||2 < \\U n — S n \\2 + \\S n — ^45 n ||2 < C. (28) 

Exploiting ((2TJ) for W 2 = W n and W\ = S n gives 

•M(W n - S n ) = M{W n ) - M(S n ) - (W n - S n , S n - A 2 S n ) 

< C{W n ) + \N(S n )\ + (||W n ||oo + ||5 n ||oo)||5 n - ^ 2 5 n ||! < C, 

and with (|28p we obtain 

\\W n - S n g = M(W n - S n ) + \\U n - AS n f 2 < C. (29) 
Now we are able to show Woo G C. From (|29p we infer that 

\\E K W n - W sh || 2 < \\E K W n - E K S n \\ 2 + \\E K S n - E K W sh \\ 2 <C+\\E K S n - E K W sh \\ 2 , 
and with S n ((p) — > W s h(ip) as n — > oo for all y> 7^ J.,-, we find 

\\EkWoo - W sh || 2 < c. 

Passing to the limit -KT — > 00 now gives W G H, and W G C follows because Woo was defined as 
weak* limit in W 1,0 °. 

5tep 3. There remains to show that Woo minimizes C. From (|2ip we infer that 

M{W n ) = M{E K W n ) + M(W n - E K W n ) + (W n - E K W n , E K W n - A 2 {E K W n )) (30) 

and 

\M(E K W n ) - M(EkWoo)\ < M(E K W n - EkWoo) 

+ \\E K W n - EkWMEkWoo - .4 2 (^Woo)|| 2 

hold for all K and n G N U {00}. Combining (|31|) with (j26|) gives 

M{E K W n ) HS^WEkWoo). 

Moreover, since E K W n - A 2 (E K W n ) is supported in \—K - 1, K + 1] for all n, we also have 

(W n - E K W n , E K W n - A 2 {E K W n )) (Woo - EkWoo, EkW^ - A 2 (E K W 00 )), 
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and passing to the limit n — > oo in (|3Q|) provides 

liminf N{W n ) > ^(EkWoo) + (Woo - E K Woo, E k Woo - A 2 {E K W OQ )), (32) 

n— >oo 

where we used that M(W n — E^Wn) > according to Lemma El On the other hand, evaluating 
(|3Up for n = oo gives 

= ^(EkWoo) + M{Woo - EkWoo) + (W^ - E K Woo, E K Woo - A 2 {E K Woo)) 

and due to M(Woo — EkWqo) — > as K — > oo we find 

WEkWoo) + (Woo - E K Woo, E K Woo - A 2 (E K Woo)) M(Woo). (33) 



The combination of (13211 and ()33p reveals 

liminf ^/"(W„) > M(Woo), (34) 

n— >oo 

and Fatou's Lemma provides 

liminf P(PF n ) >P(T^oo) (35) 

due to ^ > and since W n converges to Woo pointwise. Adding ([34^1 and (j35|) we conclude that 
Woo is in fact a minimizer of £, and Lemma [TT] guarantees that Woo solves the front equation 

(USD- □ 

We conclude with some remarks. 

1. Theorem [1] follows by combining Lemma [21 Lemma HI Lemma [TT1 Remark 1121 and Theorem 
M 

2. The assertions of Theorem [16] can be sharpened as follows. For each minimizing sequence 
(W n ) n we have equality signs in both ([M]) and ([35]) . and hence 

lim lim M(W n - E K W n ) = 0. 

K— »oo n— >oo 

This implies that there is only one interval Ji = Jj t and in turn that W n — Woo strongly 
in L 2 . In this sense each minimizing sequence obeys exactly one transition from negative 
phase to positive phase. 

3. If is increasing in [—1, +1] we can improve the existence result for fronts as follows. We 
choose r = 1 in (I) and consider the set 

C = {W € C : W' > 0}. 

Then C is an invariant set for the gradient flow of L and again one can show that C restricted to 
C attains it minimum (a proof tailored to monotone profiles is given in [HR09J). In particular, 
in this case there exist action minimizing fronts with monotone profile W. 
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4 Approximation of fronts 

In this section we illustrate the analytical results from £j3] by some numerical simulations. To this 
end we discretize the Euler scheme for the gradient flow of C, see ()24|) . as follows: 

1. Fix a finite interval [-L, +L] and introduce equidistant grid points by ipk = —L + 2kL/D, 
where k = 0...D and D € N is large. 

2. Approximate each profile W G C by the discrete vector W% = W(<pi) and impose the boundary 
conditions Wj = — 1 and Wi = +1 for i < and i > D, respectively. 

3. Replace the integrals in the definition of A by Riemann sums with respect to the <^j's. 

4. Choose A sufficiently small and initialize the iteration (|24|) with shock initial data Wi = sgn^?,,. 



graph of 4' Profile W for A = 0.9, step = 20 Residual for A = 0.9, step = 20 




Figure 3: Three examples with admissible potential as in Assumption [5] where $>' is plotted in the 
invariant interval. During the iteration the profiles W converge to a front. 

graph of 4> graph of <t graph of 4> 




Figure 4: Potentials with front parabolas for the examples from Figure[3](top row) and[5] (bottom 
row). 
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graph of 4' Profile W for A = 0.5, step = 1 Profile W for A = 0.5, step = 20 




graph of *' Profile W for A = 0.5, step = 1 Profile W for A = 0.5, step = 40 




Figure 5: Three counterexamples for the existence of action minimizing fronts. The first two 
examples satisfy the macroscopic constraints but violate the graph condition (G), so the iteration 
minimizes the action via an extending plateau. The third example violates the constraint $(+1) = 
1), so the profiles converge to one of the asymptotic states. 



In numerical simulations, the resulting iteration scheme has good convergence properties and de- 
creases the action provided that A is sufficiently small and L is sufficiently large. Three examples 
with normalized front data and $ as in Assumption [5] are shown in Figure [3l where <&' is always 
plotted over the invariant interval [-T, T). For plots of $ see Figure |H In the first example <&' is 
increasing in [—1, +1] and the front profile W turns out to be monotone. We refer to the remark 
at the end of $3] for an explanation, and to |HR09| for more examples. The other two examples in 
Figure [3] illustrate that the front profiles for non-convex are in general non-monotone. 

In order to illustrate the necessity of the graph condition (G) we present in Figure [5] simulations 
for potentials that violate this condition. In the first example the interval [—1, 1] is invariant under 
<J>' but ^ is negative in this interval. After some initial iterations the profiles exhibit an extending 
plateau at —1 < w* < 1, where w* is the minimizer of in [—1, 1] and satisfies w* = &'(w*). 
The onset of the extending plateau is a direct consequence of the energy landscape of C, see 
Remark [T2j The second simulation provides another counterexample for the existence of action 
minimizing fronts. Here \I/ has the correct sign close to ±1 but attains again a negative minimum in 
— 1 < w* < 1. As before, the profiles minimize their action by converging to the 'global minimizer' 
W = io*. Finally, the third example illustrates what happens if the asymptotic states violate the 
macroscopic constraints. More precisely, here we violate (jlOp due to 3>(+l) > $(—1), and observe 
that the profiles converge to the global minimizer W = +1. 
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